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O VI and Multicomponent H I Absorption Associated with a 
Galaxy Group in the Direction of PG0953+415: Physical 
Conditions and Baryonic Contentj^Q 

Todd M. Tripp^ and Blair D. Savage^ 
ABSTRACT 

We report the discovery of an O VI absorption system at z^^^s = 0.14232 
in a high resolution FUV spectrum of PG 0953+415 obtained with the Space 
Telescope Imaging Spectrograph (STIS). Both lines of the O VI AA 1032, 1038 
doublet and multicomponent H I Lya absorption are detected, but the N V 
doublet and the strong lines of C II and Si III are not apparent. We examine 
the ionization mechanism of the O VI absorber and find that while theoretical 
considerations favor coUisional ionization, it is difficult to observationally rule 
out photoionization. If the absorber is coUisionally ionized, it may not be in 
equilibrium due to the rapid cooling of gas in the appropriate temperature 
range. Non-equilibrium coUisionally ionized models are shown to be consistent 
with the observations. A WIYN survey of galaxy redshifts near the sight line 
has revealed a galaxy at a projected distance of 395 kpc separated by ~130 km 
s~^ from this absorber, and three additional galaxies are found within ^ 130 km 
s^^ of this redshift with projected separations ranging from 1.0 Mpc to 3.0 Mpc. 
All of these galaxies are luminous (0.6 - 4.0 L*), and two of them show the 
[O II] A3727 emission line indicative of active star formation. The galaxies with 
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[O II] emission are probably normal spirals. Combining the STIS observations of 
PG0953+415 with previous high signal-to-noise observations of H1821-I-643 with 
the Goddard High Resolution Spectrograph (GHRS), we find two O VI systems 
with Wi > 60 mA and z < 0.3 over a total redshift path Az of only 0.10. Both 
of these QSOs were originally observed to study the low z Lja lines and should 
not be biased in favor of O VI detection. If these sight lines are representative, 
they imply a large number of O VI absorbers per unit redshift, dN/dz ~ 20. 
The corresponding value of dN/dz for low z Lya hues with Wr > 50 mA is 
102 ± 16. We use this sample to obtain a first estimate of the cosmological mass 
density of the O VI systems ai z ^ 0. If further observations confirm the large 
dN/dz derived for the O VI systems, then these absorbers trace a significant 
reservoir of baryonic matter at low redshift. 

Subject headings: cosmology: observations — galaxies: halos — intergalactic 
medium — quasars: absorption lines — quasars: individual (PG0953-I-415) 

1. Introduction 

Hydrodynamic simulations of cosmological structure growth predict that when the 
initial density perturbations collapse, gas should be shock-heated to temperatures of 10^ - 
10^ °K (Ostriker & Gen 1996; Dave et al. 1999; Gen & Ostriker 1999a). The fraction of 
the gas which has been heated to these temperatures increases with decreasing redshift, 
and at the present epoch, the model of Gen & Ostriker (1999a) predicts that 47 % of the 
baryons (by mass) arc in this shock-heated phase, hereafter referred to as warm/hot gas (to 
distinguish it from the hotter gas in rich galaxy chisters which arc readily detected X-ray 
sources). This warm/hot gas prediction has not been adequately tested by observations 
because the soft X-rays emitted by gas at these temperatures are difficult to detect with 
current instrumentation, especially at lower temperatures where corrections for foreground 
absorption and emission are complicated. However, it may be possible to detect gas in the 
lower half of this temperature range via absorption fines of species such as O VI, Ne VIII, 
or Mg X in the spectrum of a background QSO (Verner, Tytler, & Barthel 1994). It is 
important to search for this warm/hot gas as part of the census of matter in the universe 
and because it could affect the formation and evolution of galaxies and galaxy groups and 
clusters (e.g., Blanton et al. 2000). 

There are some indications that warm/hot gas is present in some galaxy groups. For 
example, Mulchaey et al. (1996) have noted that ROSAT observations show that poor 
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galaxy groups which are rich in elhptical galaxies tend to exhibit X-ray emission (E > 0.5 
keV) while spiral- rich groups do not. They suggest that spiral-rich groups contain cooler 

4 X 10^ K) intragroup gas which is not easily detected in X-rays, and they predict 
that the intragroup medium of spiral-rich groups will produce absorption lines of O VI, 
but not C IV or N V because their column densities are too low at T > 5 x 10^ °K. 
Interestingly, Savage, Tripp, & Lu (1998) have recently identified a QSO absorber in the 
spectrum of H1821+643 (2;qso = 0.297) which fits this description: O VI absorption lines 
with two nearby spiral galaxies but no accompanying C IV or N V lines. However, the 
absorption could be due to the halo of the closer spiral galaxy (which is at a projected 
distance of 105 hj^ kpc) rather than the intragroup medium, and Savage et al. show 
that the O VI absorption could plausibly arise in very low density photoionized gas. The 
O VI doublet has been identified in several other intervening absorption systems^ both at 
moderate redshifts (e.g., Bergeron et al. 1994; Buries & Tytler 1996; Jannuzi et al. 1998; 
Lopez et al. 1999; Churchill & Charlton 1999) and at high redshifts (Kirkman & Tytler 
1997,1999). Composite spectra and statistical techniques have also been used to show 
that O VI absorption is present at high z (Lu & Savage 1993; Dave et al. 1998). In most 
cases it has been difficult to pin down the ionization mechanism definitively, partly due to 
the low resolution of the observations made with first-generation Hubble Space Telescope 
(HST) spectrographs, but in several cases there is evidence that the O VI systems occur 
in multiphase absorbing media. An important step in this approach to the search for 
warm/hot gas is to determine whether the O VI absorbers trace coUisionally ionized gas or 
photoionized gas. 

As part of the program described by Tripp, Lu, & Savage (1998) to study low z Lya 
absorption line systems, we have recently observed the low redshift QSO PG0953-I-415 
(^Qso = 0.239) with the E140M echelle mode of the Space Telescope Imaging Spectrograph 
(STIS). This high resolution FUV spectrum has revealed another highly ionized O VI 
absorber associated with a group of spiral galaxies, and in this paper we present our analysis 
of this particular absorbing system. In §2 we review the observations and data reductions 
including measurements of galaxy redshifts with the WIYN telescope. We present in §3 
the absorption line measurements. We constrain the temperature of the gas and examine 
its ionization in §4, and we discuss the implications of the observations in §5. Throughout 
this paper we assume Hq = 75/i75 km s~^ Mpc~^ and go = 0.0. Also, all wavelengths and 



^Strong O VI absorption lines have also been detected in "associated" absorption line systems with Zabs 
~ ^QSO (e.g., Papovich et al. 2000 and references therein). These are a rather different class of absorber 
which are often known to be very close to the QSO (Hamann & Ferland 1999). In this paper, we have 
focused our analysis and discussion on the intervening systems. 
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redshifts reported here are heliocentric, but in this direction hehocentric and LSR velocities 
are nearly identical.^ 

2. Observations 

2.1. STIS Spectroscopy of PG 0953+415 

PG 0953+415 was observed with STIS on 1998 December 4 and 1998 December 11 
resulting in a total integration time of 24478 seconds.[] The observations were obtained with 
the medium resolution FUV echelle mode (E140M) with the 0.2x0.2" aperture. Kimble 
et al. (1998) report that this STIS mode provides a resolution of R = A/AA = 46000 
(FWHM ^ 7 km s"^) with the 0.2x0.06" slit. With the 0.2x0.2" slit, the FWHM is nearly 
identical but the broad wings of the line spread function are stronger (see Figure 13.87 in 
the Cycle 9 STIS Instrument Handbook). Photoevents were recorded by the FUV-MAMA 
detector in accumulation mode with on-board orbital Doppler compensation, and individual 
exposures ranged from 2215 to 2880 seconds in duration. Short exposures of a wavelength 
calibration lamp were obtained with each individual observation of the QSO. 

The data were reduced with the software developed by the STIS Instrument Definition 
Team (IDT). Pixel-to-pixel sensitivity variations in the individual exposures were first 
corrected with a postlaunch fiatfield, and then the individual images were aligned by 
cross-correlationQ and then coadded. To prolong the lifetime of the MAMA detector, the 
position of the dispersed spectrum on the detector is changed every few weeks. This was 
done between the first observation of PG 0953+415 and the second, which is beneficial 
because when the data are aligned and coadded, any residual fixed-pattern noise not 
adequately removed by the fiatfield is reduced. This also enables the identification of 



^Assuming the Sun is moving in the direction / — 53°, b — 25° at 16.5 km s ^ (Mihalas & Binney 1981), 
^^LSR = Whciio - 0.1 km s"^ 

■^The QSO was observed for 10743 seconds on December 4 {HST archive ID no. 04X002010) and for 
13735 seconds on December 11 (i/ST archive ID nos. 04X001010 and 04X001020). 

®The images can be ahgned by cross-correlating the target images or the corresponding wavelength 
calibration lamp images. In principle direct cross-correlation of the target images is preferred because 
drift of the target in the aperture could cause shifts derived from the comparison lamp exposures to differ 
from the shifts appropriate for the target exposures. However, in this case the target images must contain 
enough sharp features to produce a sharp peak in the cross-correlation. We tried both methods, and we 
found that for this QSO observation, slightly better resolution was always obtained by cross-correlating the 
wavelength lamp exposures to determine the image shifts. 



- 5 - 



spurious features by comparison of the wavelength-calibrated spectra obtained on the two 
occasions. After coaddition of the individual images, the scattered hght correction described 
by Bowers and Lindler (1999) was applied. Then the spectra from each order were extracted 
with a simple unweighted slit with the standard height of 11 "lores" pixelsj^and smoothing 
of the background region also followed standard procedures for the STIS echelle modes. 
During the extraction, obviously hot pixels were fixed by interpolation between the adjacent 
pixels in the dispersion direction. The wavelength scale was computed with a postlaunch 
dispersion relation (which is updated when the position of the spectrum on the detector 
is changed) and put on a heliocentric basis. Finally, the data were flux calibrated and 
overlapping regions of adjacent orders were coadded with weighting based on signal-to-noise 
(S/N). The wavelength coverage is ~1 150- 1730 A with 5 small gaps between orders at A > 
1630 A. 

Samples of the final spectrum are shown in Figure |l|. We identify the lines at 1178.8 
and 1185.3 A as the O VI AA1032,1038 doublet at Zabs = 0.14232, and the corresponding 
H I Lya line is well-detected at 1388.7 A. A marginal absorption line is apparent at the 
expected wavelength of Ly(3 as well, but it is recorded in a region where the spectrum is 
quite noisy. Note that the Lya profile shows several components, the strongest of which 
is well-aligned with the O VI lines. There are no convincing alternative identifications of 
these lines. They are unlikely to be due to the ISM of the Milky Way since there are no 
resonance absorption lines at these wavelengths (see Morton 1991). The lines at 1178.8 and 
1185.3 A cannot be Ly/3,7, or 5 lines due to extragalactic absorption systems at different 
redshifts because the corresponding Lya lines are not detected at the expected wavelengths. 
Similarly, we have searched for alternative metal line identifications at other redshifts, and 
we find no convincing candidates. 

Since the O VI doublet was recorded in a spectral region which is relatively noisy, 
and since there are slightly hot spurious pixels apparent in this portion of the spectrum 
(see Figure |I|), we inspected the spectra obtained on December 4 and December 11 
independently to confirm that the O VI lines are real. Both O VI lines are detected in 
the individual December 4 and December 11 spectra, and the relative line depths and 
wavelength difference of the lines are consistent with the O VI identification. However, 
we noticed that the O VI A1037.6 line is broader in the December 4 spectrum than the 
December 11 spectrum. For the reader's inspection. Figure § compares the O VI A1037.6 
profiles derived from the 1998 December 4 and 1998 December 11 data and from all of the 



^The MAMAs are designed to support half- pixel centroiding (see Timothy 1994), and these half-pixels 
are usually referred to as "hires" pixels. This capability was not employed here since it does not enhance 
the data for this particular observation. 
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Fig. 1. — Portions of the STIS spectrum of PG 0953+415 obtained with the intermediate 
resolution FUV echelle mode (E140M). The O VI AA1032,1038 lines and the H I Lya hne 
at -2abs = 0.14232 are marked. Unrelated Milky Way lines in the selected regions are also 
indicated. In this figure, the data have been binned two pixels into one for display purposes 
only (i.e., all measurements in the tables and text were made using the full resolution 
unbinned data). 
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Fig. 2. — Comparison of the O VI A1037.6 profile at 2;abs = 0.14232 in the data obtained 
on 1998 December 11 (upper spectrum), on 1998 December 4 (middle spectrum), and in all 
of the data combined (lower spectrum). The spectra are plotted versus restframe velocity 
where v = at z^hs = 0.14232 and are offset for clarity. As in Figure 1, the data have been 
binned two pixels into one. 
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data combined. We suspect that the 1037.6 A hne is broader in the December 4 spectrum 
due to noise, perhaps due to fixed pattern noise which was not adequately removed by the 
fiatfield available when the data were reduced. Of course, this impacts the final spectrum as 
well, and the 1037.6 A line is broader than the 1031.9 A line in the final spectrum shown in 
Figure |l], though the magnitude of the discrepancy is reduced. It is alternatively possible 
that the O VI A1031.9 line is artificially narrow due to noise or an unrecognized warm pixel 
which fills in the profile somewhat - we cannot rule out or favor any of these possibilities 
with the current STIS data. Consequently, we provide below results from measurement of 
the 1031.9 A line only as well as results from measurement of both O VI lines. 

2.2. WIYN Galaxy Redshift Measurements 

As part of our continuing study of the relationship between low z Lya absorbers and 
galaxies (see Tripp, Lu, & Savage 1998), we have also measured the redshifts of galaxies 
in the ~ 1° field centered on PG 0953+415 with the fiber-fed multiobject spectrograph 
(Hydra, see Barden & Armandroff 1995) on the WIYN telescope. The observations and 
measurement techniques are described in Tripp et al. (1998), and the new redshifts will be 
provided in a subsequent paper. Here we note that this survey has revealed four galaxies 
within ~130 km s~^ of the O VI absorber at Zabs = 0.14232. The measured redshifts of 
these galaxies are summarized in Table |l] along with some of their properties including 
projected distance from the sight line (p), velocity displacement from the absorber, and O 
and E magnitudes from the revised APS digitization of the POSS I plates (Pennington 
et al. 1993). We also list Bj magnitudes calculated using the relation between Oaps and 
Bj derived by Odewahn & Aldering (1995) and the corresponding absolute magnitudes in 
Table 0. These are all luminous galaxies ranging from 0.6L* to 4.0L* (assuming M^* = 
-19.5 from Loveday et al. 1992). The velocity dispersion of this group is very uncertain 
since we do not know the redshift of the group center of mass and because of the small 
sample. However, if we assume that the O VI absorption arises at the center of mass of 
the group, then the radial velocity dispersion is roughly 100 km s~^, consistent with typical 
values observed in poor groups (e.g., Zabludoff & Mulchaey 1998). 

The spectrum of the galaxy in Table |l| with the smallest projected separation from the 
sight line (p = 395 kpc) shows strong Balmer lines in absorption as well as the other usual 
absorption lines, e.g., Ca II H & K, and well detected [O II] and H/3 emission lines with 
observed equivalent widths of 10.5±2 A and 4.5±0.5 A, respectively. However, the [O III] 
emission lines are not apparent (Ha is redshifted beyond the red end of the spectrum). 
Based on its luminosity and the Kennicutt (1992a,b) atlas of integrated galaxy spectra. 
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Table 1. Galaxies within ~130 km s of the O VI Absorber 



Galaxy 








Decl.'^ 


Oaps'^ 


(0 - E)aps' 






Redshift^ 


(Mpc) 


(km s^"*^) 


(J2000) 


(J2000) 










0.14280 


0.395 


126 


9 56 38.90 


41 16 46.6 


18.4 


1.3 


17.9 


-21.0 


0.14282 


1.00 


131 


9 57 21.41 


41 20 17.3 


19.0 


1.4 


18.5 


-20.4 


0.14257 


2.38 


66 


9 58 15.07 


41 23 23.8 


20.4 


2.5 


20.0 


-18.9 


0.14274 


3.00 


110 


9 55 14.17 


41 3 23.3 


19.1 


0.6 


18.6 


-20.3 



^Internal redshift errors are estimated to be ~50 km s ^. 

'^Projected distance to the sight line (impact parameter) assuming Ho = 75 km s"^ Mpc~^ 
and go = 0.0. The QSO coordinates (J2000) are R.A. = 9''56'"52!4, decl. = +41'il5'22'.'0. 

""Av = c(zgai - -Sabs)/ (1 + z^e^n) where 2;mean IS the mean of Zahs and Zgai. 

"^Units of right ascension are hours, minutes, and seconds, and units of declination are 
degrees, arcminutes, and arcseconds. 

"^POSS I O and E magnitudes from the revised APS catalog (Pennington et al. 1993), 
see |ittp:/ /aps. umn.edu . 

magnitude calculated from the relation derived by Odewahn & Aldering (1995): 
Oaps - 18 = 0.96(Sj - 18) + 0.52. 

^Absolute magnitude calculated using the interstellar extinction correction based on 
E{B — V) from Lockman & Savage 1995 and the i^'-correction K = 2.5 log {1 + z). 
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this implies that the galaxy is a normal Sb-Sc spiral or a peculiar SO-Sa galaxy, without 
an active nucleus. The galaxy at p = 3.0 Mpc is also an emission line galaxy with [O II], 
[O III], H/3, and even H7 seen in emission. This object has strong emission lines; the 
observed [O II] equivalent width is 37±10 A. There are no indications of an active nucleus 
(e.g., broad H/3 or large [O III]/H/9 ratios), so given its high luminosity we classify this 
galaxy as a starbursting spiral or a normal Sc-Sd spiral. The other galaxies in Table |l| do 
not have emission lines and cannot be unambiguously classified. 

It is important to note that the galaxy redshift survey we have been able to carry 
out to date is not very complete except at relatively bright limiting magnitudes, and it is 
quite possible that there are additional galaxies at 2; ~ 0.142 which are closer to the sight 
line than those in Table |I|. Consequently, we cannot comment on whether or not a dwarf 
galaxy could give rise to the O VI system, for example. As we shall discuss in §5.1, it is 
difficult to discriminate between the hypothesis that the O VI absorption originates in an 
intragroup medium and the hypothesis that it occurs in the ISM of an individual galaxy 
given the limited information currently available. More galaxy redshift measurements and 
deep imaging would be valuable. 



3. Absorption Line Measurements 

The absorption line profiles of the O VI and H I Lya lines at z^hs = 0.14232 are 
plotted versus restframe velocity in Figure ^. This figure also shows the spectral regions 
of lines of interest which are not detected. Restframe equivalent widths (W^r) of the O VI 
and H I absorption lines were measured using the software of Sembach & Savage (1992), 
but adapted to use the statistical uncertainties directly calculated from the signal and 
background counts recorded by STIS. This software accounts for errors due to uncertainty 
in the height and curvature of the continuum as well as a 2% uncertainty in the fiux zero 
point in the overall uncertainty in W,- [see the appendix in Sembach & Savage (1992) for 
details]. These equivalent widths are listed in Table |^. This software was also used to set 
upper limits on the N V, Si III, Si II, and C II equivalent widths, which are also listed in 
Table |. 

To measure column densities, we have used two methods: the apparent column density 
technique (e.g.. Savage & Sembach 1991) and Voigt profile fitting. For Voigt profile fitting, 
we have used the software of Fitzpatrick & Spitzer (1997) with the line spread functions 
for the E140M mode with the 0.2x0.2" aperture shown in Figure 13.87 of the Cycle 9 
STIS Handbook. This software provides Doppler parameters (b), velocities (f ), and column 
densities (A^) for a number of components which is subjectively specified by the user. The 
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Fig. 3. — Absorption line profiles of the O VI and H I Lja lines at ^abs = 0.14232, plotted 
versus restframe velocity. The dotted lines show the continua adopted for equivalent width 
and apparent column density measurements as well as profile fitting (see text §3). Spectral 
regions of lines of interest which are not detected are shown in the lower panels. In this 
figure, the data have been binned two pixels into one (for display purposes only). 
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Table 2. Equivalent Widths and Integrated Column Densities of the O VI Absorber at z 

= 0.14232 



Species 






log N^^ 








(A) 


(mA) 




(km s~^) 


(km s~^) 


HI 


1215.67 


423±28 


14.00±0.03 


-175 


330 


VI 


1031.93 


84±15 


14.0410:1° 


-40 


40 


VI 


1037.62 


97±15 


14.35iO:;« 


-40 


40 


N V 


1238.80 


<36^ 


<13.2'i 


-40 


40 


Si IV 


1393.76 


<m'^ 


<12.9<^ 


-40 


40 


Si III 


1206.50 


<41'i 


<12.3^ 


-40 


40 


Si II 


1304.37^ 


<27^ 


<13.3<i 


-40 


40 


CII 


1334.53 


<42'i 


<13.3'i 


-40 


40 


CIV 


1548.20 


<280f 


<13.8s 







^Restframe vacuum wavelength from Morton (1991). Oscillator 
strengths used for these measurements were also obtained from 
Morton (1991). 

''Rcstframc equivalent width integrated from w_ to w+. This 
velocity range includes all readily apparent components (see Figure 
2). 

'^Apparent column density (see text) integrated from v- to 
"^4(7 upper Umit. 

''The strongest Si 11 line at 1260.42 A is recorded in a region of the 
spectrum affected by the FUV MAMA repeller wire, so we elected to 
use the Si II A1304.37 line instead. In the wavelength range covered 
by the STIS spectrum, the Si II lines at 1190.42 and 1193.29 A are 
expected to be stronger than the 1304.37 A line, but these lines 
have unreliable oscillator strengths (see §3.2.1 in Tripp, Lu, & Savage 
1996). 

^4(7 upper limit from the Faint Object Spectrograph observations 

reported by Jannuzi ct al. (1998). 

SThe 4(7 upper limit for log A^(C IV) obtained assuming 
unsaturated absorption. 
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Table 3. Velocities, Doppler Parameters, and Column Densities of Individual Components 

of the z = 0.14232 Absorber 



Species 


(A) 


(km s~^) 


b 

(km s~^) 


log N 


HI 


1215.67 


-142±3 


121^5 


12.74±0.10 






1±2 


31±7 


13.59±0.03 






81±2 


30±9 


13.43±0.04 






163±7 


~20^ 


12.6±0.21 






207±3 


191:1' 


13.08±0.08 






269±4 


17+26 
-■-'-10 


12.81±0.12 


VI 


1031.92<= 


-3±3 


19tl' 


14.03±0.09 



OVI doublet* -2±2 221^" 14.16±0.06 



^Individual component parameters determined from 
profile fitting with the software of Fitzpatrick & Spitzer 
(1997). 

^Restframe vacuum wavelength from Morton (1991). 
Oscillator strengths used for these measurements were 
also obtained from Morton (1991). 

^Velocities are in the restframe of the O VI absorber 
with ti = km s'^ at z^^s = 0.14232. 

^This parameter is poorly constrained by the profile 
fitting software; the nominal value at the minimum is 
listed. 

^Component parameters based on fitting of the O VI 
A1031.92 line only. For the reasons noted in §2.1, this fit 
is preferred over the fit to both lines of the O VI doublet. 

^Component parameters based on fitting of both lines 
of the O VI AA1031.92, 1037.62 doublet. 
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profile fitting results are summarized in Table 0. Since the O VI A1037.6 A line may be 
artificially broadened by a noise feature (see §2.1), we present two alternative fits in Table 
^: a fit to the O VI A 1031. 9 line only and a fit to both of the O VI lines. 

In the apparent column density approach, the apparent column density per unit 
velocity, Ni,{v), is calculated directly from the apparent optical depth, 

N^iv) = Mne^)if\)-'r,iv) = 3.768 x 10'\f\)-Hn[Uv)/Iiv)] (1) 

in atoms cm~^ (km s~^)~^, where / is the oscillator strength, A is the wavelength in A, 
I{v) is the observed line intensity and /c(^^) is the estimated continuum intensity at velocity 
V, and the other symbols have their usual meanings. If the lines are fully resolved or 
are not affected by saturation, then the total column density can be obtained by simple 
integration, N = J Ng,{v)dv. We used the software of Sembach & Savage (1992) to measure 
the apparent column densities, so again, uncertainties due to the continuum placement and 
a 2% uncertainty in the fiux zero point are included in the final error bars. 

The H I Lya and O VI A1031.9 A^a(^) profiles are plotted in Figure ^, and their 
integrated apparent column densities are provided in Table 0. This table also lists 4a 
upper limits on column densities of several species of interest obtained by integrating iVa(f ) 
over the velocity range of the undetected lines. We also include a 4a upper limit on C IV 
A1548.2 from the Faint Object Spectrograph (FOS) observations of Jannuzi et al. (1998). 
Figure ^ shows that there appears to be a small velocity difference between the centroid of 
the strongest component of the H I profile and that of the O VI line. However, the H I and 
O VI velocities from profile fitting agree within their la uncertainties. More importantly. 
Figure ^ shows that the main component of the H I profile is broader than the O VI 1031.9 
A line, which suggests that the line widths are (at least partially) controlled by thermal 
motions (see §4.1). However, it is alternatively possible that the two strongest components 
of the H I profile include contributions from narrow features which are not detected in the 
O VI line. Such additional narrow lines, if present, are not well-constrained by the current 
data. 



4. Physical Conditions 

We now turn to the physical conditions of the absorbing gas. As discussed in §1, we 
are especially interested in testing the warm/hot gas prediction of cosmological simulations 
of structure growth. For this purpose we seek to constrain the gas temperature in §4.1, and 
we examine the ionization mechanism in 54.2. 
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Fig. 4. — Apparent column density profiles (see text §3) of (a) the H I Lja and (b) O VI 
A1031.9 lines plotted versus restframe velocity where w = at Zabs = 0.14232. There appears 
to be a small difference between the velocity centroids of the O VI and the main component 
of the H I profile. However, the component velocities from profile fitting agree within the la 
uncertainties (see Table ^). These Ns,{v) profiles are constructed from the full resolution, 
i.e., unbinned, data. 
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4.1. Gas Temperature 

If the line broadening is entirely due to thermal motions, then the gas temperature can 
be determined from the Doppler parameter of a given line, 

T = mh'^/2k = A(6/0.129)^ (2) 

where m is the mass and A is the atomic mass number of the element, the Doppler 
parameter b is in km s~^, and T is in °K. However, it is usually true that additional factors 
such as gas turbulence or multiple unresolved components contribute to the width of a line, 
and consequently the gas temperature from equation (2) must be treated as an upper limit. 
If lines of two or more elements with adequately different masses are available, then b can 
be expressed as 

b' = bl, + (0.129)2TM, (3) 

which can be solved for T and bnt-, the component of the broadening due to non-thermal 
motions, assuming the different absorption lines arise in the same gas and that the 
non-thermal motions (turbulence) can be assumed to have a Gaussian profile. 

The temperature upper limits implied by equation (2) and the Doppler parameters 
in Table | are T(0 VI) < 3.5 x 10^ °K and T(H I) < 5.8 x 10^ °K for the strongest 
component of the H I profile, which is also the closest component to the O VI line in velocity 
space (see Figure |). Since the velocity centroids of the O VI and the H I are quite close 
and, in fact, in agreement within the la error bars, it is plausible that these O VI and H I 
lines originate in the same gas. Adopting this assumption, we derive from equation (3) T = 
3.9 X 10^ °K and bnt = 18 km s~\ 

While such low temperatures clearly favor photoionization (see below), we must bear 
in mind that there are considerable uncertainties in the Doppler parameters. For example, 
for the strongest component of the H I profile, increasing the 6-value in Table ^ by just 
la increases the upper limit on T(H I) to 8.7 x 10^ °K, much closer to the temperature 
range where O VI can be produced by collisional ionization. This issue is exacerbated by 
the non-uniqueness of the fitted profile when dealing with multiple blended components, 
i.e., assuming a different number of components or a different mix of broad and narrow 
components can lead to substantially different 6-values. We provide an example of this in 
§4.2.2. It is also possible that the absorption arises in non-equilibrium collisionally ionized 
gas in which cooling is faster than recombination (see §4.2.2). 
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4.2. Ionization 

Given the measurements in the previous sections, can we favor photoionization or 
collisional ionization in this O VI absorber? In the Milky Way ISM, O VI likely traces 
coUisionally ionized hot gas (Jenkins 1978a,b), at least in the disk and lower halo. However, 
in extragalactic O VI absorbers, the EUV ionizing radiation field may be substantially 
harder, and the absorption might arise in very low density gas with a long path length. 
Both of these factors make photoionization more viable in the extragalactic case. We first 
consider the plausibility of photoionized models (§4.2.1), and then we test the collisional 
scenario (§4.2.2). 

4-2.1. Photoionization 

We have explored the photoionized scenario using CLOUDY (version 90.04; Ferland et 
al. 1998) with the standard assumptions, in particular that there has been time for thermal 
and ionization equilibrium to prevail (see below). The absorber is treated as a constant 
density plane-parallel slab photoionized by the extragalactic radiation from QSOs and 
AGNs, as calculated by Haardt & Madau (1996) for z ~ 0.12. We set the mean intensity at 
the H I Lyman limit to J,y(LL) = 1 x 10~^^ ergs s~^ cm~^ Hz~^ sr~^, a value in agreement 
with observational constraints (e.g., Kulkarni & Fall 1993; Maloney 1993; Tumlinson et 
al. 1999) and theoretical expectations (ShuU et al. 1999). With these assumptions, 
we varied the metallicityQ and ionization parameter {U = n^/nn = H ionizing photon 
density/total hydrogen number density, neutral + ionized) to search for models which are 
consistent with the constraints set above: the H I and O VI column densities, limits on 
various column density ratios which can be derived using the upper limits in Table 0, and 
the gas temperature. The relative heavy element abundances (e.g., [N/0]) were initially set 
to the solar values from Grevesse & Anders (1989) and Grevesse & Noels (1993), then we 
considered an alternative N/0 abundance (see below). 

The basic results of the photoionization modeling are encapsulated in Figure ^. 
This figure shows various column density ratios involving O VI and other high ions of 
interest, plotted as function of the ionization parameter (bottom axis) and riH (top axis). 
A metallicity of Z = 1/10 Zq was used to construct this model, but at low metallicities 



^°We express the linear abundance of element X relative to element Y as (X/Y) and the logarithmic 
abundance in the usual fashion, [X/Y] — log (X/Y) - log (X/Y)©, and we indicate the overall metallicity 
with the variable Z. 
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Fig. 5. — Photoionization model calculations of logarithmic ratios of the column densities of 
O VI to N V (open circles), O VI to C IV (filled squares), and O VI to Ne Vlll (open stars) 
plotted versus the ionization parameter (bottom axis) and the hydrogen number density (top 
axis). The temperature of the model as a function of the ionization parameter is also shown 
with a dotted line using the linear scale on the right axis (which is in units of 10^ °K). The 
observational lower limit on N{0 VI)/A^(N V) requires log U > -1.26, i.e. only the region 
of the plot to the right of the vertical dashed line is in agreement with this observational 
constraint. 
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the column densities scale nearly directly with Z!^ The gas temperature of the model as 
a function of U is also indicated by a dotted line in Figure |^ with the linear temperature 
scale on the right axis. Of the various species that we are able to place limits on using the 
STIS E140M spectrum (see Table only N V usefully constrains the photoionization 
model assuming [N/0] = 0; when the lower limit on A^(0 VI)/A^(N V) is satisfied, the 
other species covered in the STIS E140M spectrum are predicted to have undetectably 
small column densities. Note that the STIS spectrum does not cover the C IV doublet at 
this redshift, and we discuss the constraint set by the FOS upper limit on C IV below. 
The O VI/N V column density ratio requires log U > -1.26. Given this limit on U, the 
photoionized model requires n-n < 10^^'^ cm^^ and T > 31,200 °K. At this ionization 
parameter, the O VI ionization fraction /(O VI) = O Vl/Ototai = 0.164, and therefore 
the path length L through the constant density absorbing region must be greater than 
~460 kpc to produce the observed O VI column density, log A^(0 VI) = 14.04, in 1/10 
solar metallicity gas. With such a long path length, expansion of the universe should 
lead to an O VI line which is broader than the 1031.9 A line shown in Figure ^, and the 
photoionized model predicts a substantially larger H I column density than observed. While 
the narrowness of the O VI line may be due to noise, the incorrect A^(H I) rules out the 
1/10 solar metallicity photoionized model. To match the observed A^(0 VI) and A^(H I) 
with log U = -1.26 and J;/(LL) = 10~^^ ergs~^ s~^ cm~^ Hz~^ sr~^, the photoionization 
model requires Z ~ 1.6Zq and L ~ 25 kpc, which also alleviates the problem with Hubble 
broadening and the width of the line. While this is an uncomfortably high metallicity to 
require, recent ASCA observations of X-ray bright galaxy groups indicate that in some 
groups, the intragroup gas metallicity is ~l/3-l/2 solar (Davis, Mulchaey, & Mushotzky 
1999; Hwang et al. 1999), so this may be marginally plausible. However, in such intragroup 
gas coUisional ionization should be important. 

Are these reasonable physical properties? Observations have established that some 
QSO absorbers can be quite extended. For example, the characteristic Lya absorber "size" 
inferred from close QSO pairs is several hundred kpc at moderate redshifts (e.g., Dinshaw 
et al. 1995; Fang et al. 1996), and this is consistent with cosmological simulations 
of structure formation (e.g.. Ranch, Haehnelt, & Steinmetz 1997; Dave et al. 1999). 
However, as various authors have argued on different grounds, these large absorbers are 
likely filamentary structures composed of smaller clouds rather than a single monolithic gas 
cloud (e.g.. Ranch, Weymann, & Morris 1996; Cen & Simcoe 1997). Cen & Simcoe report 
that individual clouds have sizes of ~50 kpc. Therefore the path length required by the 



^-'^As Z approaches the solar value, the increased cooling provided by the metals affects T and, in turn, 
the ionization balance of the gas. 
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photoionized model, L ~ 25 kpc, is reasonable. The low densities inferred for the O VI 
absorber at 2;abs = 0.14232 are also consistent with cosmological simulations: according to 
equation (2) in Dave et al. (1998), the strongest component of the 2;abs = 0.14232 absorber 
is expected to have nn ~ 10~^'^ cm~^ based on its measured H I column density. Such low 
densities are also possible in the outer halo of a galaxy like the Milky Way. Murali (2000) 
has recently argued that the existence of the Magellanic Stream requires rin ^ 10~^ cm~^ 
in the Milky Way halo at a Galactocentric distance of ~50 kpc. Evidently, the size and 
density of this O VI absorber inferred from the photoionized model are plausible. 

The most problematic property of the photoionized model is the required metallicity, 
Z ~ I.GZq. Most environments in which high metallicities are expected should involve 
coUisional ionization or else should have substantially higher H I column densities. If we use 
expression (7) from Dave et al. (1999) to estimate the approximate overdensity associated 
with the 

^abs — 0.14232 absorber based on its observed A'^(H I), we find that Ph/ph ~ 
20. Based on this overdensity, the cosmological simulations of Cen & Ostriker (1999b) 
predict that the metallicity of this absorber should be Z ~ O.IZq. Similarly, to produce the 
absorption in a high metallicity region of the ISM of a galaxy like the Milky Way, the H I 
column density would have to be substantially higher than the observed A^(H I). 

We note that the H"*" recombination timescale, i:rec(H''') = l/a{T)nc where a{T) is the 
recombination coefficient, exceeds the age of the universe by a factor of several at such low 
densities, and one may wonder if modeling the gas with a photoionization equilibrium code 
is valid. However, the timescale for the gas to approach equilibrium, tgq = [^kin + ^^rec]"^ 
where tion is the photoionization timescale (see Appendix A in Dove & Shull 1994), is 
vastly shorter than the recombination timescale in the conditions considered here, and 
photoionization equilibrium is a good approximation for the low density, highly ionized gas 
(see, e.g., Vedel, Hellsten, & Sommer-Larsen 1994). 

A caveat in the photoionization models above is that we have assumed the solar [N/0] 
abundance. If instead we assume [N/0] ~ -1.5 based on the trend of [N/0] vs. metallicity 
observed in giant extragalactic H II regions (Vila-Costas & Edmunds 1993), then the lower 
limit on A^(0 VI)/iV(N V) is satisfied at a lower ionization parameter: the model requires 
log U > -1.90. However, at lower values of U the model predicts that Si III and Si IV are 
detectable, and in fact in this case the lower limit on A^(0 VI)/A^(Si III) from Table |^ 
provides a more stringent constraint requiring log U > -1.69 and < 10~^'^^ cm~^. 
Similarly, the FOS upper limit on A^(C IV) can be used to set tighter constraints on the 
photoionized model in this case: the A^(0 VI)/A^(C IV) limit requires log U > -1.33 and 
nn < 10~^-^^ cm~^. 



- 21 - 



4-2.2. Collisional Ionization 

We next consider whether the gas could be colhsionally ionized. In §4.1 we derived 
T < 5.8 X 10^ °K for the gas giving rise to the main component of the H I profile 
(the component at w ~ km s~^). According to the collisional ionization equilibrium 
calculations of Sutherland & Dopita (1993), no O VI is produced by collisional ionization at 
this temperature, and in fact O II is the dominant ionization stage of oxygen. If we assume 
that the O VI and H I absorption lines arise in the same gas with the same turbulent 
broadening, we obtain T = 3.9 x 10'^ °K, which makes collisional ionization seem rather 
unlikely. However, the absorber could be out of ionization and thermal equilibrium or it 
could be a multiphase medium, possibilities which are discussed below. 

Furthermore, as noted in §4.1, there are substantial uncertainties in the 6- values, and 
the H I Doppler parameter is consistent with gas having T ^ 1.6 x 10^ °K at the 3cr 
level. From the O VI Doppler parameter we obtain T < 3.5 x 10^ °K. These temperatures 
are more in line with colhsionally ionized gas. Furthermore, it is entirely possible that the 
H I profile contains a broad component at f ~ km s~^ superposed on a more narrow 
component. To demonstrate this, we show in Figure ^ two independent fits to the H I 
Lya line. The model profile indicated with a solid line is the initial best fit which has the 
profile parameters summarized in Table The model profile plotted with a dotted line 
used a very similar initial guess at the component parameters but was forced to include 
an additional component with the same velocity as the O VI line and b = 76 km s~^, the 
Doppler parameter of H I in gas with T = 3.5 x 10^ °K and no turbulent broadening. 
The two profile fits are nearly indistinguishable — by adjusting the parameters of the other 
components, the profile fitting code is able to compensate for the presence of the broad 
component to produce a very similar final result. The final column density required by the 
profile fitting software for the broad hot component is A^(H I) = 1.4 x 10^^ cm~^. Unless 
the spectrum has very good S/N, it would be easy to miss such a broad hot H I component. 
Therefore we cannot rule out the possibility that the gas is colhsionally ionized on the basis 
of the H I Lja profile properties. 

Figure |] is not shown in order to claim that such a broad component is present, but 
rather to show that it is allowed by the data. However, this exercise clearly shows that if 
the O VI absorption occurs in hot gas with T 3 x 10^ °K, then the absorber must be a 
multiphase medium because cooler gas is required to account for the rest of the H I Lya 
absorption including much of the absorption centered at f ~ km s~^. Given sufficient 
sensitivity, metal lines such as C III A977.02 and Si III A1206.50 should be detectable 
in these cooler phases, but this may require substantial improvement over the current 
signal-to-noise due to the low H I column densities. 
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Fig. 6. — Two independent fits to the H 1 Lya profile at 2;abs = 0.14232 (observed spectrum is 
plotted with a histogram). The solid line shows the initial best-fit with the line parameters 
in Table |^. The dotted line shows an alternative fit which started with a very similar 
initial guess at the line parameters but which was forced to include a broad component 
with a velocity and 6-value consistent with hot gas associated with the O VI absorber. This 
additional broad component is shown with a dashed line offset by 0.3 units above the data. 
It is very difficult to distinguish between the two fits, and we cannot eliminate the possibility 
that the H I profile contains a broad component due to hot gas. Note that the STIS data 
have been plotted at full (i.e., unbinned) resolution in this figure. 
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We next examine the high ion column density ratios predicted for colhsionally ionized 
hot gas. Several useful high ion column density ratios from the coUisional ionization 
equilibrium calculations of Sutherland & Dopita (1993) are plotted as a function of gas 
temperature in Figure 0. From this figure we see that if the gas is in collisional ionization 
equilibrium, then T must be greater than 2.3 x 10^ °K to satisfy the lower limit on 
A^(0 VI)/iV(N V), assuming [N/0] = 0. If [N/0] = -1.5, then T > 1.8 x 10^ °K. These are 
plausible temperatures as discussed above. 

However, the temperature range shown in Figure ^ is near the peak of the cooling 
function, and gas at these temperatures is likely to cool faster than it can recombine, 
even if it has low metallicity (Sutherland & Dopita 1993), leading to an overionized 
condition. Consequently, it may be more appropriate to compare the observed ratios to 
non-equilibrium calculations such as cooling fountain and cluster flow models (e.g., Edgar 
& Chevalier 1986; Benjamin & Shapiro 2000). Table § compares the observed lower 
limits on the O Vl/Si IV, O VI/C IV, and O VI/N V column density ratios in the z^hs = 
0.14232 absorber to the ratios predicted by non-equilibrium colhsionally ionized models 
including the cooling fountain calculation of Benjamin & Shapiro (2000), the cluster cooling 
flow model of Edgar & Chevalier (1986), the turbulent mixing layer model of Slavin, 
Shull, & Begelman (1993), and the magnetized thermal conduction interface modeled by 
Borkowski, Balbus, & Fristrom (1990). From this table we see that the cooling fountain and 
conductive interface models are fully consistent with the current observational constraints 
on this absorber given the right choice of model parameters. The turbulent mixing layer is 
marginally inconsistent with the O VI/C IV limit, but this should be confirmed with higher 
sensitivity observations of A^(C IV). However, we note that if the turbulent mixing layer 
post-mixing gas temperature is increased, then the model ironically becomes inconsistent 
with the O Vl/Si IV limit as well, presumably due to increased production of Si IV by 
self-photoionization. 

Given the available observational data, it is difficult to definitively constrain the 
ionization mechanism. It would be interesting to search for the Ne VIII A A 770.41, 780.32 
doublet in O VI absorbers. From Figures ^ and ^ one can see that in the equilibrium 
ionization models, O VI/Ne VIII ^ 1 unless the gas has rather low density (in the 
photoionized case) or is relatively hot (if colhsionally ionized). On the other hand, A^(0 VI) 
~ A^(Ne VIII) in the non-equilibrium cooling fountain model (Table Therefore detection 
of Ne VIII would provide evidence in favor of non-equilibrium collionsional ionization. 
Given sufficient z, the Ne VIII doublet will be redshifted into the bandpass of the Far 
Ultraviolet Spectroscopic Explorer or even the STIS A range. Unfortunately, ^abs = 0.14232 
is inadequate, and the Ne VIII doublet is unobservable in this system. It is possible to look 
for the S VI AA 933.38, 944.52 doublet at ^abs = 0.14232 with FUSE. However, in most 
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Fig. 7. — Logarithmic high ion column density ratios expected in hot colhsionally ionized 
gas according to the equihbrium calculations of Sutherland & Dopita (1993): O VI/N V 
(open circles), O VI/C IV (filled squares), and O VI/Ne VIII (open stars). The ratios are 
plotted versus gas temperature. The observational lower limit on N{0 VI)/iV(N V) requires 
log T > 5.36, i.e. only the region of the plot to the right of the vertical dashed line is in 
agreement with this observational constraint. 
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Table 4. Comparison of High Ion Ratios to Non-Equilibrium Collisionally Ionized Models 



Column Density 


Observed 


Cooling 


Cluster 


Mixing 


Conductive 


Ratio 


Limit^ 


Fountain*^ 


Cooling Flow'^ 


Layer'^ 


Interface" 


Vl/Si IV 


>13.8 


6.2 - 49 


66 - 180 


17 - 19 


39 - 210 


VI/C IV 


>1.6f 


2.0 - 10 


5.7 - 14 


0.7 - 0.7 


1.4-4.5 


VI/NV 


>5.5 


13 - 23 


14 - 21 


7.6 - 7.6 


4.1 - 7. 4 


VI/S VI« 




120 - 140 


92 - 104 




44 - 81 


VI/NeVIIIh 




1.4 - 2.4 









^Based on the data in Table ^ unless otherwise indicated. 

'^Cooling fountain model from Benjamin & Shapiro (2000), including self- 
photoionization, for gas initially heated to 10^ °K. 

'^Cluster cooling flow model from Edgar & Chevalier (1986) for the isochoric and 
Ttr = 2.2 X 10^ °K cases. 

'^Turbulent mixing layer ratios calculated by Slavin et al. (1993) for entrainment 
velocities of 25-100 km s~^, a post-mixing gas temperature of 10^ ° °K, and no dust 
depletion (see their Table 4A). 

''Thermal conduction front model from Borkowski et al. (1990) with an age of 10^'^ 
years and magnetic field inclinations ranging from to 60° relative to the front normal 
(see their Figure 6). Note that in this model, these ratios increase or remain roughly 
constant as the age exceeds 10^-^ years. 

^Based on the 4ct upper limit on the C IV equivalent width from Jannuzi et al. 
(1998), assuming the linear portion of the curve of growth applies. This limit is less 
conservative than the limits in the third column of Table |^ because it does not account 
for the uncertainty in the continuum placement. 

VI has resonance lines at 933.38 and 944.52 A which can be observed with the 
Far Ultraviolet Spectroscopic Explorer (Sembach 1999). Unfortunately, at this redshift 
one of the S VI lines may blend with the Ar I 1066.66 A resonance line. 

''The Ne VIII resonance lines are at 770.41 and 780.32 A and therefore cannot be 
observed at the redshift of the O VI absorber studied in this paper. However, FUSE 
will be able to look for this doublet in somewhat higher redshift O VI systems such 
as the one studied by Savage et al. (1998), so we include it here for purposes of 
illustration (see text). 
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models this doublet is predicted to be substantially weaker than the O VI lines (see Figures 
land 0and Table 1). 

5. Discussion 

5.1. Association with a Galaxy Group 

The fact that the z^hs = 0.14232 O VI absorption line system is associated with a group 
of galaxies strongly indicates that the absorber is an intervening system rather than an 
intrinsic absorber which was ejected or somehow accelerated to high displacement velocity 
by the QSO. While this may seem a trivial conclusion since the O VI is displaced from 
the QSO redshift by ~23,000 km s~^, during the last few years observations of absorption 
variability have established that some highly ionized intrinsic QSO absorbers are separated 
from zqso by such large velocities (e.g. Hamann, Barlow, & Junkkarinen 1997) and yet 
are relatively narrow (i.e., not traditional broad absorption line outflows). Therefore, it is 
important to find evidence that a given O VI absorber is indeed intervening even if Av is 
large. Of course, this association with a galaxy group does not necessarily indicate that 
the absorption arises in the intragroup medium; it could be due to gas within one of the 
galaxies in the group. We briefly discuss some possibilities below. 

The close proximity of the O VI system to a galaxy group provides a theoretical 
prejudice in favor of collisional ionization. The presence of galaxies requires that 
intergalactic gas has collapsed in this region of space, and simple arguments (e.g., §3 in Cen 
& Ostriker 1999a) suggest that substantial shock-heating probably occurred as a result. 
Therefore coUisionally ionized gas is expected in the vicinity (if not along the pencil-beam 
probed by the QSO). It is interesting that this absorption system and associated galaxies 
fit the prediction of Mulchaey et al. (1996) discussed in §1: a group which is possibly 
spiral-rich and has associated O VI absorption, as expected based on their postulated 
coUisionally ionized intragroup medium. However, we really are not sure that this is a 
spiral-rich group (or even a bound group) with only four known galaxies, and it is also 
possible that the O VI absorption arises in the gaseous halo of a single galaxy. However, 
models of galaxy gaseous halos (e.g.. Mo & Miralda-Escude 1996) also usually produce 
O VI absorption in a coUisionally ionized hot phase. Similarly, models which produce QSO 
absorption line gas in supernova-driven winds from dwarf galaxies (e.g. Wang 1995) also 
involve substantial shock-heating and coUisionally ionized gas. 

The kinematics of the absorption also provide useful information about the possible 
origins of the system at z^hs = 0.14232. A fundamental question is whether the absorption is 
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due to the intragroup medium or to the ISM of a single galaxy that happens to intercept the 
line-of-sight. The H I Lja profile has the "leading edge asymmetry" that various authors 
(e.g., Lanzetta & Bowen 1992) have discussed as the signature of a moderately edge-on 
rotating disk. Of course, this is not a unique interpretation of such profile asymmetry, but 
it provides some evidence in favor of the single-galaxy interpretation. The simplicity of the 
O VI profile also favors this interpretation. The velocity dispersion of poor galaxy groups is 
typically one to a few hundred km s~^ (Zabludoff & Mulchaey 1998), so one might expect 
the O VI to be spread over a larger velocity range if the absorption is due to the intragroup 
medium. Here, though, we must recognize that the observed O VI profiles are noisy and we 
only detect O VI at the velocity of the strongest component of the H I profile. If the hot 
gas is concentrated at the center of the group, then this may be the only location along 
the line-of-sight where we have sufficient sensitivity to detect it, and there may be O VI 
absorption at other velocities which has fallen below our detection threshold. To sort out 
the various possibilities, it would be very helpful to obtain additional STIS observations to 
improve the signal-to-noise and search for O VI at other velocities. 

5.2. Number Density and Cosmological Mass Density 

One means to test the warm/hot gas prediction of cosmological simulations (§1) 
is to compare the number of O VI absorbers observed per unit redshift, dN/dz, to the 
number statistically predicted from many random pencil-beams through the cosmological 
simulations. Is the observed dN/dz consistent with the number predicted by the 
cosmological models? We can also estimate the mean cosmological mass density traced by 
the O VI systems at low z. 

For these purposes, we combine the STIS observations of PG0953+415 with the GHRS 
observations of H1821+643 from Tripp et al. (1998).Q We emphasize that the manner 
in which these QSOs were selected should not bias the sample to enhance the number of 
O VI systems detected compared to sampling many random directions. Both QSOs were 
originally observed to study the relationship between weak Lya clouds and galaxie^ at 
z < 0.3. These two objects were selected for this project simply because they are among the 
brightest known QSOs with z > 0.2, criteria which were required to substantially increase 



^^After this paper was completed, Tripp, Savage, & Jenkins (2000) carried out an analysis of new STIS 
echelle observations of H1821+643 with the E140M mode, and we refer the reader to that paper for the O 
VI number density and cosmological mass density derived from an independent data set. 

^^HST program IDs 6155 and 7747, see Tripp et al. (1998) for details. 
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the sample of weak Lya absorbers with a minimal amount of HST time. No consideration 
was given to factors which might indicate an enhanced likelihood of detecting the O VI 
doublet when the targets were selected. 

To estimate dN/dz of the O VI systems, we define a sample of O VI lines with > 
60 mA for both lines of the doublet]^ and set the maximum absorber redshift, Zmax? for 
each sight line to exclude any absorbers within | Av \ < 5000 km s~^ of the QSO redshift 
and thereby avoid contaminating the sample with associated/intrinsic absorbers which are 
close to the QSO. One might argue that 5000 km s~^ is insufficient (see above). However, 
we find that the two O VI systems in the final sample are associated with galaxy groups, 
which suggests that these are indeed intervening absorbers. The lower redshift cutoff for 
each sight line, z^am, was determined by the lowest wavelength in the observed spectrum 
(with a small buffer to ensure that a line would be recognized if at that A) in the case of 
H1821+643 and by the wavelength at which the S/N is unacceptably low in the case of 
PG0953+415. These criteria resulted in a sample of two O VI absorbers within a total 
redshift path of Az = 0.100 (after a significant correction of 0.067 for regions of the spectra 
in which we cannot detect either of the O VI linesQ because they would be blocked by 
strong ISM or extragalactic lines from other redshift systems). Therefore, dN/dz ~ 20 for 
O VI systems with > 60 mA at ^ < 0.3. Using the confidence limits from Gehrels (1986) 
for a sample of two absorbers, we derive 4 < dN/dz < 63 for these O VI absorbers at the 
90% confidence level. 

For comparison, Tripp et al. (1998) find dN/dz = 102 ± 16 for H I Lya lines with 
Wr > 50 mA at z < 0.3, and Weymann et al. (1998) report dN/dz = 33 ± 4 for Lya 
absorbers at z = with W,- > 240 mA based on a large sample at 2; < 1.5. In the case of 
low to moderate redshift Mg II absorbers, dN/dz = 0.97 ± 0.10 for > 300 mA (Steidel 
& Sargent 1992) and dN/dz = 2.65 ± 0.15 for > 20 mA (Churchill et al. 1999; see 
also §6 in Tripp, Lu, & Savage 1997). Evidently, these weak low 2 O VI systems have a 
substantially larger cross section and/or covering factor than the Mg II absorbers. Similarly, 
the stronger O VI absorbers at higher redshifts are less common: Buries & Tytler (1996) 
report dN/dz = 1.0 ± 0.6 for O VI systems with > 210 mA at < ^abs > = 0.9, and 
similar results (with smaller uncertainties) are derived from the larger sample provided by 
the FOS Quasar Absorption Line Key Project (B. Jannuzi & R. Weymann 2000, private 
communication). This comparison with the O VI dN/dz derived from FOS data must be 



^''The 3(7 detection limit throughout the region where O VI absorbers can be detected is 60 mA or better 
in the spectrum of PG0953+415 and 50 mA or better in the spectrum of H1821+643. 

^'"^Both of the O VI Hnes were required to faU in unblocked regions of the spectra so that the doublet could 
be securely identified. 



interpreted carefully, however, because the FOS samples are dominated by lines which 
are substantially stronger and at substantially higher redshifts than the O VI absorbers 
discussed in this paper. With a larger sample of low z weak O VI absorbers and smaller 
uncertainties in their number density, comparison of dN/dz to the space density of objects 
such as dwarf galaxies may provide insight on the nature of the O VI systems. 

Next we estimate the baryonic content of the O VI absorbers, expressed as the 
cosmological mass density fib(0 VI), following previous analogous calculations for damped 
Lya systems (e.g., Lanzetta et al. 1991) as well as O VI absorbers (Buries & Tytler 
1996). To estimate the density of baryons in the gaseous component of the universe traced 
by O VI, we require information about the metallicity of the gas and the O VI ionization 
fraction, /(O VI) = O Vl/Ototai- In collisional ionization equilibrium, /(O VI) peaks 
at ~0.2 (Sutherland & Dopita 1993), and similarly low peak fractions are predicted by 
non-equilibrium collisional models (e.g., Shapiro & Moore 1976; Benjamin & Shapiro 2000). 
The O VI ion fraction is not much larger at peak value in photoionized gas (see appendix). 
Therefore we will adopt /(O VI) ~ 0.2 to set a lower limit on the baryonic content of 
the O VI absorbers, and the following calculation is relatively independent of how the gas 
is ionized or whether or not it is due to the intragroup medium. For the mean cosmic 
metallicity of the O VI absorbers, it is less clear what value to adopt, but to set a lower 
limit on f2;,(0 VI) we should set (0/H) to a high but plausible value. As noted in §4.2.1, 
Davis et al. (1999) and Hwang et al. (1999) have derived metallicities of 1/3-1/2 solar 
for intragroup gas in several X-ray bright groups. Therefore we will initially use 1/2 solar 
metallicity for the calculation of ^2^(0 VI) and then discuss how it scales with (0/H). The 
mean cosmological mass density in the O VI absorbers, in units of the current critical 
density pc, can be estimated as 



where /i is the mean atomic weight (taken to be 1.3), (0/H)o vi is the assumed mean 
oxygen abundance by number in the O VI absorption systems, mn is the mass of hydrogen, 
NiiO VI) is the total O VI column density and AXj is the absorption distance interval 
(Bahcall & Peebles 1969) probed to the ith QSO, 



blocked by strong lines. Combining the PG0953-I-415 STIS data and the GHRS observations 



^^Over the redshift range probed by the sight hues to PG0953+415 and H1821+643, results are insensitive 
to the value assumed for q^. 





(5) 
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of H1821+643 from Tripp et al. (1998) with Zmin and z^^a.^ set to the same values used for 
the derivation of dN/dz, we obtain ^^(O VI) ;^ O.OOOG/ifg^ assuming the mean O abundance 
is 1/2 solar. This is a lower limit not only because /(O VI) and (0/H) were set to their 
approximate upper limits, but also because we have applied an equivalent width cutoff to 
define the sample; if O VI absorbers with IVj. < 60 mA significantly increase J2iNi{0 VI), 
then the true ^2^(0 VI) will be higher. Note that ^2^(0 VI) is inversely proportional to 
(0/H). Decreasing the metallicity to 1/10 solar, for example, increases the lower limit on 
the O VI absorber baryon content to Jlfe(0 VI) ^ O.OQSh^^. 

To demonstrate the level of uncertainty in the cosmological mass density estimate due 
to small number statistics, we can recalculate ^^(O VI) using an alternative expression 
analogous to equation (9) from Rao & Turnshek (2000), 



Pec/(o VI) vh;ovi v^^^y + 

where < A^o vi > is the mean O VI column density of the absorption systems in the 
sample and again we have assumed go = 0. The advantage of this alternative expression for 
flb{0 VI) is that we can employ the Gehrels (1986) small sample statistics to estimate the 
uncertainty in dN/dz and then propagate the uncertainty into the estimate of Qh{0 VI). 
With dN/dz = 20l^| and the other parameters set to the values assumed above, we find for 
the 1/10 solar metallicity case ^5(0 VI) ;^ 0.003lo'oo2 ^75^ (error bars are la uncertainties). 

Bearing in mind that there is still considerable uncertainty in the lower limit on 
f2fe(0 VI) due to the small sample, small rcdshift path, and uncertain mean metallicity, 
this preliminary estimate suggests that the O VI absorbers may indeed harbor a significant 
fraction of the baryons in the universe at low z. The lower limit assuming (0/H) = 1/10 
solar is comparable to the cosmological mass density of stars, H I, and X-ray emitting galaxy 
group and cluster gas at low redshift (Fukugita, Hogan, & Peebles 1998), for example. 



6. Summciry 

The paper is summarized as follows. 

(1) We have observed the low redshift QSO PG0953+415 with the E140M mode of 
STIS, and this has revealed an O VI absorption line system at z^hs = 0.14232. This O VI 
absorber is highly ionized: multicomponent H I absorption is also detected at this redshift, 
but no other species are detected in the STIS spectrum including the N V doublet. Si III, 
and C II. 
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(2) We have measured galaxy redshifts in the field of the QSO using the WIYN 
telescope, and there are at least four galaxies within ~130 km s~^ of the O VI absorber 
with projected distances ranging from 395 kpc to 3.0 Mpc. Two of these galaxies appear to 
be spiral galaxies. 

(3) A review of the observational constraints shows that we cannot definitively assert 
that the gas is collisionally ionized or photoionized, but photoionization requires an 
uncomfortably high metallicity which is inconsistent with theoretical expectations. If the 
gas is collisionally ionized, it is likely that it is not in equilibrium or that it is a multiphase 
absorber. Non-equilibrium coUisional ionization models are consistent with the observations. 

(4) Combining the STIS data on PG0953+415 with the high S/N low resolution GHRS 
observations of H1821+643 from Tripp et al. (1998), we identify two intervening O VI 
systems over a redshift path I\z of only 0.10. This implies that dN/dz for O VI systems 
with W-c > 60 mA and z < 0.3 is ~20 with a large uncertainty due to the small number of 
systems so far detected. This represents a large density for metal line systems since dN/dz 
for LycK absorbers with > 50 mA and z < 0.3 is 102 ± 16 (Tripp et al. 1998). We stress 
that the sample should not be biased in favor of O VI detection. If further observations 
confirm that dN/dz is as large as 20 for O VI systems, then these absorbers may be an 
important baryon reservoir at low redshift, although this depends on the metallicity of the 
gas. If the mean metallicity is 1/2 solar, then Qb{0 VI) ^ O.OOOGhj^. However, if the mean 
metallicity is 0.1 solar, then ^^5(0 VI) ^ O.OOShj^, which is comparable to the baryonic 
content of other known constituents of the low z universe such as galaxies and gas in galaxy 
clusters. 
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digital format, and Buell Jannuzi and Ray Weymann for private communication regarding 
the number density of O VI systems in the FOS Quasar Absorption Line Key Project. 
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A. Maximum Ionization Fraction of O VI 

In §5.2 we adopted /(O VI) = 0.2 to set a lower limit on the baryonic content of the 
O VI absorbers. Intuitively, /(O VI) = 1.0 would seem to set a more conservative lower limit 
on fih(0 VI). However, as noted in §5.2, O VI is not a preferred ionization stage of oxygen, 
and in collisionally ionized gas /(O VI) peaks at ~0.2 in equilibrium and non-equilibrium 
calculations. 

To evaluate the maximum O VI ionization fraction which can be attained if the O VI is 
created by photoionization, we have calculated with CLOUDY an O VI ionization fraction 
grid as a function of the temperature and ionization parameter of a parcel of gas assuming 
the Haardt & Madau (1996) UV background radiation aX z = 0.12. To show how /(O VI) 
depends on T and U , we operated CLOUDY in a different mode than that employed in 
§4.2.1: instead of allowing CLOUDY to determine the gas temperature given the various 
heating and cooling processes which are important (as was done in §4.2.1), we fixed the 
temperature and ionization parameter at particular values in each grid cell and calculated 
the ensuing /(O VI) for a grid with 3.0 < log T < 7.0 and -3.5 < log U < 0.5. The resulting 
grid is shown in Figure § 

As expected based on previous calculations, when the ionization parameter is low 
enough so that photoionization is unimportant and O VI only has a significant ionization 
fraction in gas which is hot enough to be collisionally ionized. Figure || shows that /(O VI) 
has a maximum value of ~0.2 at log T = 5.5. On the other hand, when the ionization 
parameter is high and photoionization dominates. Figure indicates that /(O VI) could 
be somewhat larger if the gas is very cool. However, this is rather unlikely to occur because 
the photoionization by the UV background which creates the O VI will also heat the gas to 
log T ;^ 4. For example, the CLOUDY model in Figure ^ which satisfies the O VI and 
N V constraints requires log T ^ 4.5 (see §4.2.1), and with this temperature constraint 
the maximum /(O VI) is 0.28. Therefore /(O VI) = 0.2 is a reasonable value for placing 
a lower limit on the baryonic content of O VI absorbers regardless of whether the gas is 
photoionized or collisionally ionized. 



tl 



ion 




Fig. 8. — Dependency of the O VI ionization fraction, /(O VI) = O Vl/Ototai, of a parcel of 

gas on the ionization parameter U and gas temperature T. assuming the gas is photoionized 
by the UV background calculated by Haardt & Madau (1996) at 2; ~ 0.12 (see §4.2.1). It 
is important to note that U and T were held at fixed values for the calculation of /(O VI) 
in each cell of this grid to illustrate how the O VI ionization fraction varies with these 
parameters. However, the low temperature region of this plot is not physically realistic 
because the photoionization which ionizes oxygen to this degree will also heat the gas to log 
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